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ABSTRACT 


This  thesis  investigates  the  performance  degradation  resulting  from 
multitone  interference  on  orthogonal,  noncoherent,  frequency-hopped,  M-ary 
frequency-shift  keyed  receivers  (FH/MFSK)  where  the  effect  of  thermal  and  other 
wideband  noise  is  not  neglected.  The  multiple,  equal  power  jamming  tones  are 
assumed  to  correspond  to  some  or  all  of  the  possible  FH  M-ary  orthogonal 
signaling  tones.  Furthermore,  the  channel  is  modeled  as  a  Ricean  fading  channel, 
a  possibility  precluded  when  thermal  noise  is  neglected.  Both  the  signaling  tones 
and  the  multiple  interference  tones  are  assumed  to  be  affected  by  channel  fading. 
Both  band  and  independent  band  multitone  interference  are  considered. 
Performance  is  evaluated  by  obtaining  a  union  bound  on  the  probability  of  bit 
error,  and  receiver  performance  is  compared  with  exact  results  for  band  multitone 
interference  of  a  noncoherent  FH/MFSK  receiver  under  comparable 
circumstances.  Except  for  the  case  of  Rayleigh  fading  of  the  signal,  the  union 
bound  is  very  tight  for  those  cases  that  can  be  compared  with  exact  results.  The 
advantages  of  the  union  bound  approach  are  twofold.  First,  the  union  bound 
approach  yields  a  solution  that  is  far  less  computationally  intensive  than  that 
obtained  with  the  exact  approach.  Second,  the  union  bound  approach  allows 
numerical  results  to  be  obtained  for  interference  conditions  that  are  not  amenable 
to  exact  analysis,  such  as  independent  multitone  interference  of  FH/MFSK. 
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I.  INTRODUCTION 


A.  SCOPE  OF  THESIS 

The  goal  of  this  thesis  is  to  determine  the  electronic  counter-counter  measures 
(ECCM)  potential  of  a  non-coherent  frequency-hopped  M-ary  frequency-shift  keying 
(FH/MFSK)  communication  system  under  conditions  of  worst  case  hostile  electronic 
counter  measures  (ECM)  and  fading  channels.  The  performance  of  the  FH/MFSK 
receiver  will  be  investigated  for  frequency-nonselective,  slowly  fading  Ricean  channels 
and  hostile  multitone  interference. 

The  receiver  considered  uses  a  quadratic  detector  in  each  M-ary  branch  of  the 
receiver  for  signal  detection.  Both  band  and  independent  multitone  interference  are 
considered  and  both  jammer  and  signal  are  assumed  to  be  affected  by  the  fading  channel, 
where  it  is  also  assumed  that  the  channel  fading  need  not  necessarily  affect  the  signaling 
tones  and  the  interference  tones  in  the  same  way.  The  cases  of  Rayleigh,  Ricean,  and  no 
fading  are  all  considered. 

Performance  is  evaluated  by  obtaining  a  union  bound  on  the  probability  of  bit 
error.  The  results  are  compared  with  the  exact  results,  examined  in  [1],  to  determine  the 
accuracy  of  the  union  bound  approach.  Then  the  union  bound  is  used  to  obtain  probability 
of  bit  error  for  situations  where  the  exact  approach  is  either  very  difficult  or  even 
impossible  to  obtain. 


B.  IMPORTANCE  OF  A  SPREAD  SPECTRUM  COMMUNICATION  SYSTEM 

The  benefit  from  investigating  the  performance  of  a  spread  spectrum 
communication  system  lies  in  the  extent  to  which  spread  spectrum  communications  are 
currently  used.  Their  application  to  any  military  communication  system  aids  in  achieving 
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secure,  low  probability  of  intercept,  and/or  low  probability  of  detection  information 
exchange  between  friendly  units  as  well  as  reducing  the  ability  of  hostile  elements  to  jam 
friendly  communications.  Other  applications  of  spread  spectrum  are: 

1 )  The  use  of  either  frequency-hopping  or  direct  sequence  spread  spectrum  in 
multiple  access  communications  (as  in  cellular  phones  where  vehicles  communicate  with 
a  central  station)  in  which  a  number  of  independent  users  are  required  to  share  a  common 
channel  without  an  external  synchronization  mechanism. 

2)  The  application  of  spread  spectrum  techniques  in  wireless  local  area  networks 
(LANs).  Wireless  LANs  are  very  important  for  the  establishment  of  a  network  for  mobile 
forces  that  circulates  among  the  cooperating  units  various  kinds  of  information,  such  as 
text  messages,  video,  or  audio.  Wireless  LANs  are  also  very  popular  in  commercial 
communications  and  are  used  for  wireless  communication  of  laptop  computers  with  a 
control  station. 

In  any  case,  the  potential  of  frequency-hopping  spread  spectrum  to  reject 
interference  of  either  unintended  or  intended  transmissions  with  the  communicating 
signal  is  of  crucial  importance.  The  problem  considered  in  this  thesis  is  the  evaluation  of 
the  performance  of  a  frequency-hopped  system  in  order  to  determine  the  system's 
capabilities  and  limitations. 

In  chapters  2,  3,  and  4,  the  union  bound  approach  to  performance  analysis,  a  brief 
description  of  frequency-hopping,  and  the  notation  used  in  the  thesis  are  discussed.  In 
chapters  5  and  6,  receiver  performance  with  independent  and  band  multitone  jamming, 
respectively,  are  analyzed.  In  chapter  7  numerical  results  for  both  independent  and  band 
multitone  jamming  are  presented. 


II.  IMPORTANCE  OF  UNION  BOUND  APPROACH 

The  union  bound  approach  simplifies  the  mathematical  development  by  deriving 
an  expression  for  the  upper  bound  on  the  probability  of  symbol  error  that  is  always  much 
simpler  than  the  exact  form.  There  are  usually  two  approaches  used  to  evaluate  the 
probability  of  symbol  error  of  a  communications  system. 

In  the  first  approach,  known  as  the  indirect  approach,  the  probability  of  symbol 
error  is  obtained  using  the  relation 

Ps=l-Ps  (1) 

where  Ps  is  the  probability  of  not  making  a  symbol  error.  Noting  by  X,,X2,...XM  the 
random  variables  representing  the  output  of  each  branch  of  the  receiver  for  a  signal  with 
M  symbols,  we  have: 

Ps=l-Ps=l-Pr(X1>X2nX1>X3n...nX1>XM|l).  (2) 

The  calculation  of  this  probability  is  difficult,  especially  in  the  case  where  we  deal  not 
only  with  gaussian  noise  but  also  with  jammer  tones  affecting  the  signal.  This  will 
become  particularly  obvious  later  on  when  we  analyze  receiver  performance  with 
independent  multitone  interference.  In  that  case,  the  joint  probability  distribution 
functions  required  to  solve  the  problem  using  the  exact  solution  approach  make  it 
impossible  to  obtained  a  closed  form  expression  for  the  probability  of  symbol  error  for  a 
system  using  more  than  two  symbols. 

However,  if  we  use  the  union  bound,  or  direct  approach,  the  solution  works  as 
follows,  providing  an  upper  bound  to  the  probability  of  error  calculation: 


Ps  =Pr(X,  <X2uX,  <X3u...uX,  <XM|1)  <  Pr(X,  <X2|l)  +  Pr(X,  <X3|l) 
+...+Pr(X1<XM|l) 


We  see  that  in  place  of  joint  probability  distribution  functions  (pdfs)  we  now  have  the 
sum  of  individual  pdfs.  In  fact,  the  individual  pdfs  can  be  grouped  in  groups  of  identical 
pdfs  occurring  with  a  certain  multiplicity,  so  that  the  solution  is  a  linear  combination  of  a 
few  pdfs. 

Summarizing,  the  union  bound  approach  provides  a  simple,  fast  method  to 
investigate  communication  and  jammer  signal  scenarios  that  are  very  difficult,  if  not 
impossible,  to  explore  with  the  exact  result  approach  as  described  in  [1].  By  using  the 
union  bound,  we  can  increase  our  knowledge  of  communication  systems  performance  for 
a  larger  variety  of  jamming  conditions.  However,  since  this  evaluation  is  a  bound,  its 
accuracy,  or  tightness,  needs  to  be  tested  so  that  we  can  comment  on  its  effectiveness  and 
limitations  as  a  tool.  In  fact,  in  the  following  pages  it  will  be  shown  that  in  some  cases  the 
union  bound  approach  is  so  accurate  that  it  provides  virtually  identical  results  with  the 
exact  solution,  while  in  some  other  cases  the  results  obtained  display  an  anomalous 
behavior. 
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III.  BACKGROUND 

The  main  idea  of  frequency-hopping  spread  spectrum  communications  is  to 
provide  protection  against  a  hostile  jammer  by  increasing  the  bandwidth  that  the 
information  signal  occupies  much  more  than  the  minimum  required.  By  increasing  the 
frequency  range  over  which  a  system  effectively  operates,  we  force  the  jammer  to  spread 
its  power  over  a  wider  frequency  band  and  thus  be  less  effective  from  a  "per-frequency" 
signal  corruption  capability  point  of  view. 

The  main  idea  of  frequency-hopping  spread  spectrum  is  based  on  the 
multiplication  of  the  signal  prior  to  transmission  by  an  intermediate  frequency,  generated 
by  a  frequency  synthesizer,  of  the  form 

fi  =  f,  +  (i  - 1)  •  Af ft  where  i=l  ,2,...N  (4) 

where  N  is  the  maximum  number  of  possible  frequency  hop  bins,  Afft  is  the  separation 
between  the  carrier  frequencies  of  adjacent  bins,  and  T  changes  pseudo-randomly  every 
Tc  seconds.  By  so  doing,  the  entire  spectrum  of  the  signal  transmitted,  is  shifted  from  its 
carrier  frequency  fc  to  the  new  carrier  frequencies: 

fCi=fc+fr+(i-l)-Afft.  (5) 

The  result  is  that  a  new  data  modulated  carrier  that  is  shifted  from  one  frequency  to  the 
next  is  created.  From  (5)  we  see  that  there  are  N  different  frequency  hop  bins,  each  of 
bandwidth  Afft  .  The  value  of  i  (and  thus  fCj )  is  changed  periodically  according  to  some 
predetermined,  but  apparently  random  to  a  third  party  observer,  noiselike  spreading  code 
called  either  the  pseudo-random  or  the  pseudo-noise  (PN)  sequence. 


Generally,  frequency-hopping  is  used  with  noncoherent  modulation  formats  such 
as  noncoherent  MFSK.  Such  a  signal  is  described  by: 

s(t)  =  V2Accos[27i(fs+(m-l)-Af)t  +  0i]  (6) 

where  m=l,2,...M,  0;  is  the  i*  symbol  phase,  and  Af  is  the  frequency  separation  between 
each  of  the  M  signaling  tones.  Multiplying  this  signal  by  a  signal  at  the  intermediate 
frequency  f;  =  f,  +  (i  - 1)  •  Afft ,  we  get: 

s'(t)  =  2a/2Ac  cos{2Tc[fs  +  (m  -  l)Af  ]t  +  9j }  cos(27tfjt)  = 

=  V2AC  cos{27r[fi  +  fs  +  (m  -  l)Af  ]t  +  9{ }  +  V2AC  cos{27i[f;  -  fs  -  (m  - 1) Af  ]t  +  G; } 

(7) 

The  carrier  frequencies  of  the  first  part  are 

fi+fs+(m-l)Af=f1+(i-l)Affh+fs+(m-l)Af  (8) 

and  is  smallest  for  i=l  and  m=l;  i.e.,  fs  +f, .  The  carrier  frequencies  of  the  second  part 
are 

fi-fs-(m-l)Af=fi+(i-l)Affh-fs-(m-l)Af  (9) 

and  is  largest  for  i=N  and  m=l;  i.e., 

f1+(N-l)A4-fs.  (10) 


If  the  smallest  frequency  of  the  first  part  is  greater  than  the  largest  frequency  of  the 
second  part,  then 

2fs  >2f,  +2B  +  (N-l)Affh  =>fs  >f,  +B  +  (N-l)Affh  II  (11) 

where  B  is  the  required  guardband  above  and  below  the  high  and  the  low  frequency 
signaling  tones  respectively.  In  this  case,  the  signal  s'(t)  can  be  high-pass  filtered  to 
remove  the  frequency  difference  contribution,  and  the  frequency-hopped  signal  becomes: 

V2AC  cos{27i[f,  +fs+(m-l)Af]t  +  0i}  (12) 

One  obvious  advantage  is  that  a  hostile  jammer  needs  to  know  the  pseudo-noise 
sequence  of  the  communicator's  transmitter  in  order  to  jam  the  specific  hop  bin  where  the 
transmitter  operates  at  each  time  instant.  To  overcome  a  follower  jammer,  the  signal  must 
hop  to  a  new  carrier  frequency  before  the  jammer  is  able  to  complete  its  tracking  process. 

Another  advantage  is  due  to  the  fact  that  the  power  spectral  density  of  the 
frequency-hopped  signal  is  identical  to  the  that  of  the  conventional  signal  in  a  specific 
hop  bin.  However,  since  the  signal  hops  from  bin  to  bin,  and  assuming  that  the 
probability  that  any  bin  is  occupied  is  equal  to  (1/N),  the  average  power  spectral  density 
is: 

PSD=^£SFSK(^c=fc,)-  (13) 

^       i=l 

Hence,  the  signal  power  spectral  density  is  lowered  by  a  factor  of  N  on  average. 

Before  proceeding  we  need  to  distinguish  between  the  two  main  tone  jammer 
methods  used  against  communication  systems.  The  first  is  band  multitone  interference 
and  the  second  is  independent  multitone  interference.  In  each  case,  we  assume  that  the 


multiple  jamming  tones  are  all  of  equal  power.  In  band  multitone  interference,  the 
multiple  equal  power  jamming  tones  are  distributed  randomly  across  the  entire  frequency 
hop  bandwidth,  but  the  number  of  jamming  tones  placed  within  each  jammed  frequency 
hopped  band  is  specific.  In  independent  multitone  interference,  the  multiple  equal  power 
jamming  tones  are  distributed  in  a  completely  random  fashion,  meaning  that  the  number 
of  jammer  tones  within  a  frequency  hop  bin  can  vary  from  zero  to  a  maximum  equal  to 
the  total  number  of  signal  tones  corresponding  to  each  symbol  the  communicator  can  use. 
We  also  assume  that  each  jamming  tone  corresponds  exactly  to  one  of  the  possible 
M  x  N  signaling  tones.  From  the  communicator's  viewpoint,  this  is  the  worst  case. 

Before  proceeding  to  the  analysis  of  the  receiver  performance  with  multitone 
jamming,  the  notation  used  in  the  remainder  of  the  thesis  is  discussed  in  the  next  chapter. 


IV.  FH/MFSK  SYSTEM,  CHANNEL  MODEL,  AND  SYMBOLOGY  USED 

The  FH/MFSK  communication  system  examined  is  assumed  to  have  N  non- 
overlapping  FH  bins.  The  bandwidth  of  each  one  is  B  Hz,  where  B  is  the  bandwidth  of 
the  non-hopped  MFSK  signal.  So  the  total  frequency-hopped  bandwidth  is  augmented  to 
be  greater  than  or  equal  to  B-N.  We  also  express  by  M  the  modulation  order  where 

M  =  2k  (14) 

and  k  is  the  number  of  bits  per  M-ary  symbol.  Since  the  number  of  signal  tones  within  a 
hop  bin  is  M  and  the  total  number  of  bins  is  N,  the  result  is  a  total  of  M-N  possible  signal 
tone  locations. 

The  symbol  duration  Ts  is  related  to  the  hop  duration  Tc  by: 

TC=KTS  (15) 

where  K=l,2,3... .  The  symbol  and  hop  rates  are  the  reciprocal  of  the  corresponding 
durations,  so 


Rs=^  (16) 


andRc=^,  (17) 


respectively.  Symbol  rate  is  related  to  hop  rate  by 


R,=K-RC  (18) 

Using  the  frequency  hopped  signal  (12),  we  note  that  the  M  symbols  are  assumed 
to  be  equally  likely,  both  fj  and  fs  are  taken  to  be  integer  multiples  of  the  symbol  rate  Rs, 
and  the  frequency  difference  between  signal  tones  Af  is  taken  to  be  an  integer  multiple  of 
the  symbol  rate.  Hence, 


Af  =  p-Rs  (p  integer). 


(19) 


To  incorporate  the  idea  of  fading  channels,  we  model  the  amplitude  Ac  as  a  Ricean 
random  variable.  This  includes  Rayleigh  fading  as  a  specific  case. 

The  transmitter  and  receiver  are  assumed  to  function  in  perfect  synchronization, 
and  the  receiver  after  dehopping  the  signal  uses  a  bank  of  M  quadrature  detectors  where 
the  integrator  time  constants  are  normalized  to  the  symbol  duration  for  notational 
convenience  (Figure  1). 

The  existing  noise  interference  is  expressed  as  AWGN  with  a  flat  average  power 

No 
spectral  density  defined  as  — .  The  total  power  of  the  jammer  is  expressed  as  PJt  ,  and 

the  jammer  delivers  q  equal  power  interfering  tones  such  that  the  power  of  a  single 
interfering  tone  is  equal  to 


p,= 


(20) 


The  subscripts  used  in  the  description  of  the  Ricean  random  variables  are  'c'  for 
the  signal  parameters  and  'j'  for  the  jammer.  Thus,  the  pdf  of  Ac  is 


M«.) 


<T2. 


exp 


a2c+a2^ 


2a: 


Io 


^a2c-aV 
a2c     . 


U(ac) 


(21) 
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where  u(»)  is  the  unit  step  function,  I0(»)  is  the  modified  Bessel  function  of  the  first  kind 

and  order  zero,  a2c  is  the  power  of  the  direct  signal  component,  and  2a2c  is  the  power  of 
the  diffuse  signal  component  of  the  respective  tone.  The  total  average  power  of  either  a 
Ricean  faded  jammer  or  signal  tone  is  equal  to  the  mean  squared  value  of  the  amplitude  ac 
or  aj5  respectively,  and  is 

I*  =a2i+2a2i  (i  =  c,j).  (22) 

The  average  power  is  assumed  to  remain  constant  from  hop  to  hop.  It  is  convenient  to 
define  the  signal-to-noise  power  ratios  of  the  direct  signal  components  as 


Pi=4i     (i  =  c>J)  (23) 

a  n 


and  the  diffuse  signal  components  as 


2(T2i 

4i=-T^      (i-cj)  (24) 

O    n 


where  a2n  is  the  noise  power. 

For  orthogonal  signaling,  the  average  probability  of  bit  error  is  related  to  the 
average  probability  of  symbol  error  by  [2] 


p>=^.p„ 


11 


and  the  average  energy  per  bit  is  related  to  the  average  energy  per  symbol  by 
Es=log2(M)Eb=kEb.  (26) 
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V.  FH/MFSK  RECEIVER  PERFORMANCE  WITH  INDEPENDENT 

MULTITONE  INTERFERENCE 


With  independent  multitone  interference,  the  multiple,  equal  power  jamming 
tones  are  distributed  randomly  across  the  entire  frequency-hopped  bandwidth.  Two 
possible  cases  exist  concerning  the  signal  and  jammer  tones  within  a  specific  frequency- 
hop  bin.  In  the  first  case  there  is  no  interference  at  all  in  the  bin  (the  jammer  tones  are 
elsewhere),  and  the  probability  of  error  is  the  same  as  for  a  conventional  MFSK  signal  in 
the  presence  of  AWGN  only.  In  the  second  case  there  is  jammer  interference  in  the  bin. 
Since  a  bin  is  composed  of  M  different  signal  locations  and  a  potential  jammer  may 
place  one  or  more  of  its  tones  to  the  aforementioned  signal  locations,  this  case  is  further 
divided  into  two  subcases.  In  the  first  subcase  the  signal  and  jammer  tones  are  on 
different  branches,  while  in  the  second  subcase  one  jammer  tone  coincides  with  a  signal 
tone.  This  gives  a  total  of  three  individual  cases  that  are  treated  below. 

Case  1: 

In  the  following,  we  express  by  R,,  R2,...RM  the  random  variables  representing  the  outputs 
of  branches  1,2,. ..M,  respectively,  of  an  MFSK  receiver.  Without  loss  of  generality,  we 
assume  that  the  symbol  T  is  transmitted.  In  order  not  to  have  an  error,  the  output  of 
branch  one  must  be  greater  than  the  outputs  of  all  other  branches.  The  opposite  situation 
expresses  the  error  condition.  As  in  equation  (2),  the  probability  of  symbol  error  provided 
by  the  union  bound  is 

Ps  =  Pr(R,  <R2uR,<  R3u...uR,  <  RM|l) 

<Pr(R,  <R2|l)  +  Pr(R,  <  R3|l)+...+  Pr(R,  <RM|1) 
<(M-l)Pr(R,  <R2|1) 

and  the  last  equation  is  true  since  the  individual  probabilities  are  equal. 
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Case  2a: 

Suppose  that  we  have  L  jammer  tones  in  the  hop  bin  where  the  signal  is.  By  Si  we  express 
the  existence  of  a  jammer  tone,  and  the  subscript  after  j  corresponds  to  the  specific  branch 
where  the  jammer  tone  is  present.  For  example  S}   means  that  a  jammer  tone  is  present  in 

branch  one.  To  simplify  the  discussion  for  this  case,  suppose  that  we  have  two  jammer 
tones  (L=2).  Now 

Pr(error|l,Sj2,SJ3)<Pr(R,  <  R2|l,Sj2)  +  Pr(R,  <R3|l,Sjj)  +  Pr(R]  <R4|l)+...+Pr(R]  <RM|1) 
Since  Pr(R,  <  R4|l)  =  Pr(R,  <  R5|l)  =...=  Pr(R,  <  RM|1)  = 

=  Pr(error  with  no  interference  in  signal  branch)  (28) 

and  Pr(R,  <  R2 11,  SJ2 )  =  Pr(R,  <  R3 11,  SJ3 ) .  (29) 

then  (28)  reduces  to 

Pr(error|l,  SJ2 ,  SJ3 )  <  2  Pr(R,  <  R2 11,  S h  )  +  (M  -  3)  •  Pr(R,  <  R4 11)  (30) 

Generalizing  for  L  jamming  tones  per  hop  bin,  1<L<M-1,  we  get 

Pr(error|l,XSJi+i)<L-Pr(R,  <R2|l,SJ2)  +  (M-L-l)-Pr(R]  <RM|1)  (31) 

1=1 

Case  2b: 

As  with  case  2a,  we  begin  by  considering  two  jammer  tones.  In  this  case 

Pr(eiror|l,SJi,SJi)<Pr(R1  <R2|l,SJi,SJ2)  +  Pr(R]  <R3|l,SJi)  +  Pr(R1  <k4Ji,SJI)+..- 

+  Pr(R,<RM|l,SJ]) 

(32) 
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As  before,  this  can  be  simplified  to 

PrCerrorMj,  ,SJ2 )  <  Pr(R,  <  R2  lUS^  ,Sj2 )  +  (M  -  2)  •  Pr(R,  <  R2  \\,Sh ) .  (33) 

Generalizing,  for  L  jamming  tones  per  hop  bin,  1<L<M-1,  we  get 

Pr(error|l,SJ|,Xs.+i)  =  (L-l)-Pr(R1<R2|l,SJi,SJ2)  +  (M-L)Pr(R1<RM|l,SJi) 

1=2 

where  it  is  understood  that  the  summation  is  zero  for  L  =  1 . 

(34) 

Now  the  probability  of  symbol  error  with  at  least  one  jamming  tone  per  hop  bin  (L>1)  is 
given  by: 


Pr(error)=  Pr(error|l,2X+i ) .  Pr(l,£sju|  )  +  Pr(error|l,Sj|  ,£sji+i )  •  P^LS^  ,£sjm ) 

1=1  1=1  1=2  1=2 

(35) 
Since  we  have  L  jammer  tones  and  M  signaling  tones  per  jammed  frequency-hop  bin, 


PKI^JX,)^  (36) 


L  T 

and  Pr(l,YS:    )=  1 .  (37) 

So  the  probability  of  symbol  error  for  L  jamming  tones  per  hop  bin,  L>1,  is  given  by 
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P,(L)< 


M 


+   1 


(L-l)-Pt(R1<RafcSA.SJl)+(M-L).Pr(R1<R1(|l>SJl)] 

L-Pt(r,  <R2|1,S.  )  +  (M-L-l)Pr(R,  <RM|l)] 


+ 


(38) 


We  now  look  for  the  probability  of  having  L  jamming  tones  per  hop  bin.  This 
leads  to  the  expression  for  the  total  probability  of  bit  error.  For  independent  multitone 
interference,  and  since  we  assume  q  jamming  tones  randomly  interfere  with  a  total  of 
N-M  possible  signal  tones,  l<q<N-M,  the  probability  that  a  frequency-hop  bin  contains 
one  interfering  tone  is 


N-M 


(39) 


The  probability  that  a  frequency  hop  bin  contains  a  second  interfering  tone  is 


q-1 


N-M-l 


(40) 


The  probability  that  a  frequency  hop  does  not  contain  an  interfering  tone  is 


1- 


N-M 


(41) 


The  probability  that  a  frequency  hop  bin  does  not  contain  a  second  interference  tone  is 


1- 


q-1 

N-M-l 


(42) 


16 


From  (39)  -  (42),  we  derive  the  probability  that  a  frequency  hop  bin  contains  no 
interference  tones  as 


l— LUi~     q 


N-M^   V      N-M-l 


M- 


^      N-M-(M-1)J     tf\      N-M-i 


Similarly,  the  probability  that  there  is  one  jamming  tone  in  a  frequency  hop  bin  is 


M-N 

q 


I      N-M-l 


1- 


q-1 


1- 


q-1 


N-M-2J V      N-M-(M-lV 


■nu 


M-Nl=f-V      N-M-i 


(44) 


In  the  same  manner,  we  define  the  probability  that  there  are  two  jamming  tones  in  a 
frequency  hop  bin  as 


q  q-1 


M-N    N-M-l 


liVM-N-i 


1.rjt^Ti.r,_.  o-2 


N-M-l/  V      N-M-2 

q-2    < 
M-N-J" 


<  q-2 

1- 


k      N-M-(M-l) 


nu 


i  =  l 


(45) 


Combining  (43),  (44),  and  (45),  we  get  the  probability  that  L  jamming  tones  are  in  a  hop 
bin,  where  1<L<M-1,  as 
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Pr(L  jamming  tones  per  hop  bin)  = 
q-(L-l)-l 


q-1 


1- 


MN-(L-l)-l 


1- 


M-N    M-N-i 
q-(L-l)-l 


q-(L-l) 

M-N-(L-l) 


M-N-(M-l) 

L-l         „       : 


liM-N-i 


M-l 


nu- 


i=L 


q-(L-l)-r 

M-N-i     J 


(46) 


The  situation  where  L=M  must  be  treated  separately.  Looking  at  the  two 
individual  products  multiplied  together  in  (46),  we  see  that  the  first  product  has  L  terms, 
corresponding  to  the  probability  that  L  out  of  the  M  signals  of  the  bin  are  jammed,  while 
the  second  product  has  the  remaining  M-L  signals  of  the  bin  which  are  not  jammed.  Since 
when  L=M  all  the  signals  are  jammed,  for  that  case  we  need  only  the  first  product  term. 
So  the  probability  of  having  L=M  jamming  tones  in  a  hop  bin  is 


L-l 


Pr(L  =  M)  =  J  ] 


q-i 


i=0 


M-N-i 


(47) 


Now,  the  total  probability  of  symbol  error  for  0<L<M  is  given  by 


Pr(error)  =  Pr(L  =  0)  •  Ps(l)  +  Pr(L  =  l)  •  Ps(l|L  =  l)+...+  Pr(L  =  M)  •  Ps(l|L  =  M) 


(48) 


The  maximum  number  of  jamming  tones  we  can  have  in  a  bin  is  M  if  M  <  q 

and  q  if  q  <  M.  This  can  be  represented  as  min(q,M).  By  combining  equations  (46), 

(47),  and  (48),  the  probability  of  symbol  error  is  given  by  the  relation 
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M-l 


p.*  ni- 


i=0 


N-M-i 


(M-l)Pr(R,  <R2|l)  + 


min(M-l,q) 


+ 


+ 


s 


L-l 

n 


q-i 


L=1     LUM-N-i, 


M-l 


i  =  L 


w- 


q-(L-l)-!' 

M-N-i    J 


P,(L) 


rL-l 

n 

l  i=0 


q-i 


M-N-i 


7  -P,(M) 


(49) 


The  probability  Ps(M)is  obtained  by  letting  L=M  in  (38).  This  yields 

PS(M)  =  (M  - 1)  •  Pr(R,  <  R2  |l,Sj  ,Sj  ) .  Based  on  the  above  probabilistic  description,  it 

is  apparent  that  we  need  to  obtain  the  following  probabilities  for  the  MFSK  receiver: 

Pr(R,  <R2|1),  Pr(R,  <R2|l,Sj|),  Pr(R,  <R2|1,S-  ), 

andPrCR,  <  R2|l,SJiSSJ2). 

A.  Pr(R,  <R2|l) 

For  noncoherent,  orthogonal  BFSK,  Pr(R,  <  R2|l)  is  equivalent  to  the  probability  of 
bit  error: 


Pb(ac)  =  ^-e  2'2" 


1 

2 


(50) 


N 


where  a  n  =  — — 


(51) 


For  noncoherent,  orthogonal  MFSK,  the  only  thing  that  changes  is  that  Tb  is  replaced  by 
T: 
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,       N„ 

<?   n  =  "=T" 

s 


(52) 


To  account  for  the  existence  of  a  Ricean  fading  channel  we  multiply  equation  (50)  with 
(21),  the  Ricean  pdf,  and  then  integrate  over  all  possible  values  of  ac .  This  yields 


ao  oo    i 

Pb  =  JPb(ac)-fAe(ac)dac  =  j-exp 


f         ~2      \ 


{    2a2nJ   a  2 


exp 


a2c  +ac2 
2a2c 


c     c 


V   gS    J 


U(ac) 


2c, 


•exp 


-a, 


2o\ 


|acexp 


1  1 


n       a, 


^a  a    » 

c     c 

V  a2c  V 


da. 


(53) 


Let  x  =  a2c  =>  dx  =  2acdac ,  and  (53)  becomes: 


Pb=^-.exp 


a, 


2a2c 


exp 


J    o 


1  1 

+  — : 


c\ 


a. 


V  cr  c  J 


da. 


(54) 


We  then  use  the  identity  [3]: 


»      m+n 

Jx~-e-ax-Jn(2pVx)dx  = 


m!(3nexp 


a 


m+n+l 


'fill 

V  a  J 


(55) 


where  3nm(«)    is  the  Laguerre  polynomial  and  is  defined  [3]: 
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?n 


m 


'fill       y    (-1)*3   fm  +  nl  (P'l 
V  a  J      p  =  0    p!      Im-pJ 


a  ; 


(56) 


In  our  case  m=n=0  and  the  whole  polynomial  evaluates  to  unity.  In  place  of  "a"  we  have 


f    1  1 


~2    +  ~2 


ja 
while  in  place  of  "B"  we  have  — f- .  So, 


Pk   = 


1 

4a2c 

exp 

1 

~ 

4a2c 

•exp 

I 

2  +  2 

a2c 

a2c 
2?7 


a 

2a1 


j  exp(-  ax)  •  J  0  (2(3  vx  Jdx  = 

0 

2  V  0~   2,-r2      ^ 
2CT    CT    n 


2<T    cC    n 

.a2c  +a2ny 


exp 


4a, 


~      2     _L    ~2 


exp 


a2c 


2a' 


a2n 


1- 

I      a  c  +a  n J 


exp 


2  +  2- 


a2c 


2a\  +CT  n 


.(57) 


a  n 


a2n 


Using  (23)  and  (24),  we  have  the  definitions: 


a2c 


pc  =  direct  signal-to-noise-ratio  = 


2(5' 


£,c  =  diffuse  signal-to-noise-ratio  = 

<7    n 


Applying  these  definitions  to  (57),  we  get 


P^R^R^^r-exp 


2  +  ^( 


(58) 
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B.Pr(R1<R2|l,Sj1) 

This  probability  of  error  is  the  same  as  that  derived  in  section  4.  A  except  that  a 
jammer  tone  is  on  the  signal  branch.  Since  signal  and  jammer  tone  exist  on  same  branch, 
we  can  consider  that  they  constitute  a  new  kind  of  (altered)  signal  existing  on  the  branch. 
So  the  probabilistic  analysis  yields  the  same  result  with  the  only  difference  being  that  the 

"old"  signal  power  a2c  is  replaced  by 

a2c+a2j+2acaJcos0j.  (59) 

In  other  words,  the  power  of  the  signal  and  jammer  tones  add  together  but  not 
necessarily  in  phase.  We  then  integrate  with  respect  to  Gj  after  multiplying  by  its  pdf  to 

eliminate  the  dependence  on  0.  The  phase  angle  9  j  is  assumed  to  be  uniformly  distributed 

in  [0,2tt]  . 

Without  jamming,  the  mean  values  of  the  random  variables  representing  the  outputs  of 

the  in-phase  and  quadrature  components  of  the  first  branch  are  given  by: 

XM  =  V2ac  cos(G)  (60) 

Xq_  =  V2ac  sin(e)  (61) 

and  in  this  case  the  output  of  the  first  branch  is  given  by: 

X2j,  +  X2q,  =2a2c(cos2(e)  +  sin2(e))  =  2a2c  (62) 


This  results  in  the  probability  of  error: 
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Pr(R,<R2|l)  =  -exp 


'    aVO 


V     2NJ 


(63) 


With  jamming  we  have  for  the  in-phase,  quadrature  components,  and  their  squared  sum 
as  the  final  branch  output  the  following  relations,  respectively, 


XSi  =  V2ac  cos(e)  +  V2aj  cos(e;) 


(64) 


Xqi  =  V2ac  sin(e)  +  V^  sin(e;) 


(65) 


and  X2i,  +X2qi  =2 


a2c  +a2j  +2^^.00^-0'^ 


(66) 


We  define  9=  =9  —  0]  as  the  phase  of  the  jammer  with  respect  to  the  received  signal  to 
get  the  final  form  of  equation  (59).  Considering  the  equations  above,  we  can  understand 
why  we  substitute  the  previous  signal  power  a2c  in  Pr(R,  <  R2  |l,Sj1 )  by  the  composite 

power  a2c  +  a2j  +  2acaj  cosfGj  J .  Substituting  (59)  into  (63),  multiplying  by 
f9  (0  ■  j  =  — ,  and  integrating  over  9  j ,  we  get 


1     27t 
Pr(R1<Ra|l,Sjl)«~.J«p 


a2c  +a2j  +2acajcos[0Jj 


2N, 


d0; 


(67) 


We  then  use  the  identity 
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1     2n 
I0(z)  =  ^  •  jexp[±  z  •  cos((|)  -  G)]d(j) 


(68) 


to  obtain 


Pr(R,<R2|l,Sj1)  =  --exp 


a2c   < 


2c\ 


1+^ 

a  cy 


•  I 


a  c    a 


C  ~J 


°w2n    a; 


=  —  •  exp 


-^( 


2azn 


a2c+a2j) 


■I, 


C  J 


v  cj  n  y 


(69) 


However,  this  probability  of  error  is  conditional  on  both  aj  and  ac .  To  eliminate  the 

dependence  on  ac,  we  multiply  by  its  pdf  (Ricean)  and  then  integrate  for  all  possible 
values  of  ac.  After  multiplying  (69)  by  (21)  and  integrating,  we  get: 


Pr^R^Sj,^ 


•exp 


2a2n 


(a2,)   "J 


-J    o 


ic--7--exp 
o  c 


2g2„ 


(a2c) 


•exp 


a2c  +q2c 
2ct2c 


C  J 

V  a2n  J 


w^K 


2g' 


exp 


2a2n 


(a2,) 


•exp 


a 


2a2c 


ac-exp 


g    c  +CT    n 

2a  ca  n 


(a2c) 


l 


V  a2„  V 


a     a 

c         c 


V     G' 


da. 


We  now  have  the  identity  [3]: 


(70) 
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Jx-exp[-ax2]-In(px)-Jn(yx)dx  =  — -exp 


2        ..2 


4a 


>.(! 


(71) 


The  correspondence  of  the  terms  in  (70)  with  the  ones  in  the  identity  are: 


J-a< 


ct2c+g2„     0       aj 
a=  9fT2      2     »P  =  -I-»  Y=      2 

iO    C(T    n  CJ    n  G    c 


So  the  probability  of  error,  after  eliminating  the  dependence  on  ac,  is 


/  \  cr2n  f     1    a2j+a2c>l      f      ar  ^ 

Pr(R1<R2|l,SJi)=    ,   2         2Vexp---      J  -Io  /        -aj  .      (72) 


To  eliminate  the  dependence  on  a^ ,  we  multiply  by  its  pdf  (Ricean)  and  then  integrate  for 
all  possible  values  of  a^.  The  Ricean  pdf  for  aj  is  given  by: 


fAjM=;i"exp 


a2j+aj2 

2a2j 


•I, 


(a-  -di.\ 


.*? . 


u  a . 
J 


(73) 


Multiplying  equation  (72)  by  (73),  we  get: 


G2n 


2(a2c+a2„) 


I'        a' 


exp 


f    1    a2J+a2c^ 


^       2     CT    c  +0"    nJ 


-j-  exp 
^a  j 


'    a2j+a2j^ 


2a2j     ) 


°l    ~2  ,  2  "j    I       X0 

v  a  c  +a  n 


a.    -I, 


V     C72j     J 
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G2n 


2(g2c+ct2„)    ct2j 


•exp 


1    G2ja2j  +a2ja2c  +a2ca2j  +  a2jG2n  +  ct2jG2c  +  a2ja2n  ' 


(a2c  +a2n)a2j 


I, 


A       U-a^ 


0     _2 


a- 


Gc+CTn         V 


J        J 


^  cr  j  y 


a2n 


exp 


2   a2j(a2c+a2„)   a2j       ^    2     (a2c+a2n)a2j  J 


2  V\ 


1    a2j+(a2c+a2n) 


exp 


2    _2     \A 


22  /       2         2  2         2     "\ 

1    a  j(x  c  +  (a  ja  c  +oc  jcr  nJ 
2"  (a2e+a2n)a2j 


a. 


-CTC  +G'n 


aj  J  -lo 


^a     a  N 

j      j 


v   a  j   ) 


(74) 


To  integrate  (74)  with  respect  to  aj  we  again  use  identity  (71).The  correspondence  of  the 
terms  in  (74)  with  the  equivalent  ones  of  the  identity  are: 


a  =  — 


\_  (T2j+(a2c+g2n) 
2     (a2c+a2n)cr2j 


,   P  = 


a. 


G  c  +a  n 


a; 


->y=y-i 


For  the  argument  of  the  Bessei  function,  we  have: 


a. 


gy 

2a 


<J   c  +C7    n 


J-^2" 


l_  a2j+(a2c+a2n) 

2     (a2c+a2n)a2j 


=  J"3 


acaj 


2  2 

CT    j  +0"    c  +0"    n 


So  the  final  Bessei  function  has  the  form:  I, 


ac«j 


VCT2j  +CT2C  +G2nJ 


(75) 
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Now,  —  = 
2a 


\ 


r,2  ^"j 

a    J  +     ~~2 2 2 

1  va  j+a  c+a  ny 

2  (a2c+a2„)a2j 


(a2c+a2n)a2J 
a2j+a2c+a2n 


So  the  terms  appearing  in  front  of  the  integral  of  (74)  reduce  to: 


2  2  2 

G    j  +CT   c  +G    n 


•(CT2c+CT2n)-G2j 


a2n 


a2n 


2   G2j(cr2c+a2n)J     2(a2j+a2c +a2n) 


(76) 


Finally,  the  argument  of  the  exp(»)  function  that  the  integral  of  (74)  yields  is: 


P2-Y2       \<j\+o\ 


a, 


4a 


J_  g2j+(q2c+cr2n) 
2     (a2c+CT2„)a2j 


and  combining  this  with  the  leading  exponential  functions  we  get  the  final  exponential 
function  to  be  of  the  form: 


a. 


a  c  +a  n 


a 


"J 


4- 


1  a2J+(a2c+a2n) 

2  (ct2c  +a2n)a2j 


1  a2ja2c +(a2ja2c +a2j<j2n) 

2  (a2c +a2n)a2j 


1       (a2c+a2j) 
y-J+a2      '*    (77) 

Z      (J    j  +G    c  +  <^    n 


Combining  (75),  (76),  and  (77)  as  required  by  (71)  we  get: 
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Pr(R,<R2|l,SJ,)  = 


a2n 


2       ,    „.  2 


•exp 


1        a  c  +a  j 


2(a2J+a2c+a2n)       *\    2    a2j+a2c+a2n 


•I, 


(         acaj  ^ 

Va2j  +a2c  +a2ny 


Substituting 


(78) 


a2c  a2j 

Pc  —  ^.2        '     Pi  ~~   _2       '     Sc  — 

a  n  o  „ 


^j  ~    -2 


2ct2c 
a2n 


in  (78),  we  get 


Pr(R1<R2|l,SjO  =  -^.exp(--^ 


^j+^ 


^j+4c+2 


•I 


0 


2VpTp7 


(79) 


We  begin  by  designating  V,  and  V2  to  be  the  random  variables  representing  the 
outputs  of  branch  1  (where  the  signal  resides)  and  branch  2  (where  the  jammer  tone 
resides),  respectively.  Their  pdfs  are  given  by 


fV  KMe! 


2a' 


•exp 


(  «  2  ^ 

Vj  +2a  c 


n 


2cr 


n   ; 


aK 


n 


(80) 


f     (v2|l,aj)  =  — y-exp 

2  -^a   n 


(  ^  2\ 

v2+2a  J 

2a2n 


a2n 


(81) 
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Multiplying  (80)  by  (21)  and  integrating,  we  get 


\h) 


a2c 


f     «   ^ 


2        _2 


2a  n   a 


exp  -      2      •  exp 


|ac-exp  -a' 


2a 


a2n  +2a2( 


^    2a  „/ 


2    _2  *0 


2a  „a 


acV^7 


k  a  n  ; 


v  a2c  > 


da. 


We  then  use  the  identity  given  in  (71)  where 


(82) 


a  = 


a2n  +2-a: 


2a  na 


P  = 


V2vT 


2    _2  '       K  —  2        ' 

a  n 


Y  =  J-2" 
a  < 


(83) 


Now  we  have 


(3-Y 

V2Viac 
~2    ~2 

a  na  c 

2a 

a  n  +  2a  c 

a2n  +2a: 


2         2 

a  na  c 


(84) 


2v,      a2c 


P2-Y2 

a  a      a  c 

2v!a4c  +a4na2c 

4a 

2(a2n  +2a2c) 
a  na  c 

2(a2n  +2a  c)o2n02c 

(85) 


Hence,  using  (83),  (84),  (85)  in  (71),  we  can  evaluate  (82)  to  obtain 


«*(viH)- 


2(a2n  +2a2c) 


exp 


1     v,+2a2/       (  V2v7-a( 

7       \  '  M 


2    a*n  +2a"< 


2  ->      2 

a  n  +2a 


C> 


(86) 
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Similarly, 


f     I       IC     \  l  i       l       V2+2a2J 


2>  ' 

•  I 


a2n  +2ct2j 


(87) 


Using  the  random  variable  transformation  R2i  =  V, ,  we  get 


f*.(R'|l)=R~^gexp 


1  r,  +2a2c    |      [    >/2r,  -a 

2  a  n  +2a  cy 


a  n  +2a  cy 


(88) 


and  similarly, 


fR,  R2|s 


,)- 


exp 


1     r2+2aj 


2    A 


(a2n+2a2J)      *\    2   a2B+2a2j 


•L 


V2r,  -a 


2     **j 


a2n  +2a2j 


(89) 


Since  R,  and  R2  are  independent  random  variables,  the  joint  pdf  of  Rj  and  R2  is  equal  to 
the  product  of  the  individual  pdfs: 


fRiR2(r],r2H,SJ2)  =  fRi(r1|l)-fRi(r2|SJ2) 


(90) 


00   00  00 

So  Pr  R,  <R2|1,SJ2)=  JJfRi  -fR2dr2dri  =  JfRi (r, |l) ■ 

Or,  0 


JfRj(r2|Sj2)dr2  dr,        (91) 


For  the  inner  integral, 
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00  00 

ft.  M*  )*i -fer^TT  "* 

ri  ri 


r2  +2a2j 
2(a2n+2a2j) 


■L 


V^<T    n  +2CT    jy 


(92) 


Let  x  = 


x-ya2n  +2a2j 
A/a2n+2a2J  """""  A/a2n+2a2J   '       a2n+2a2J  "       CT2„+2a2j 


dx  =    ,  =  and 


(93) 


Also  to  transform  the  lower  limit  of  integration,  when  r2  =  r, ,  x  =    .        l  (94) 

>/ct2„  +2a2j 


Using  (92),  (93),  and  (94)  to  (91),  we  get 


00 

Pr(R,<R2|l.S.)=       J      x-exp 


■v/o2„+2ct2j 


2a2, 
x^     a2n  +2a2j 
~2~  2 


V2  -a, 


A 

CT2n+2a2j      J 


dx 


(95) 


Now  compare  the  integral  above  with  the  Marcum's  Q-function  [2],  [4]: 


Q(a,P)=  J 


x  •  exp 


x2+a2> 
2      > 


I0(ax)dx, 


(96) 


and  notice  that  if  we  let 


a 


a  = 


=r>P 


^a2n  +2a2j  '         f\+2^j 
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Then  (95)  is  given  by 


( 


Q(a,P)  =  Q 


V^x2" 


Va2n+2a2j  '7a2„+2a2j 


(97) 


Substituting  (97)  in  (91),  we  get 


—  •exp 


r2,+2a2c    1  ,  (    v^a, 


2(a2„+2a2c). 


2      "\ 

r  i 


■I, 


a  n  +2a  cy 


.^/o5,  +2a2j  '^/a2n  +2a2j 


(98) 


Consider  the  following  identity  [5]: 


a2    R^    R, 

—  •exp 

0       ^2      GlJ      &    1 


HZ* 


a2,+R2,' 
2a2,     . 


■M^-ldR,- 


a2, 


2  2 


If  we  let 


l-Q 


a2, 


a22 


2  2     '1/      2  2 

^(J   l+O   2      \  G    \  +G   2 


a  2 


/ 


a2i  +c22 


a22 


a2, 


.2      ,    _2      »V  _2      ,    _2 


y\G   \  +G   2     \G   i  +  CT   2y 

(99) 


o~2  =  ya2»  +2a2j  ,  a2=ai-4l,  R,  =  r,  ,  a,  =  Vtf2n  +  2a2c  ,  a,=>/2-<xc 


then  (99)  can  be  used  to  evaluate  (98)  with  the  result 
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Pr(R1<R2|l,Sj2)  = 


■2g2< 

.2       ,    . 


ct„  +2a: 


2(a2n+a2c+a2j) 


1-Q 


a 


a' 


2  2  2     '1        2  2  2 


+ 


a2n  +2a2j 


2(a2n  +a2c+CT2j) 


a 


a2c 


2  2  2     '1/       2  2  2 

CTn+CJc+CTj      yCTn+CTc+C7jy 


Using  (23)  and  (24),  we  get 


Pr(R1<R2|l,Sj2)  = 


^-■Q 


i  +  4< 


2  +  4c+^j 


i-Q 


V^pT        V2^" 


V2  +  Sc+W2  +  Sc+Sj 


+ 


2  +  Sc+Sj 


V^pT       V^; 


^+47^'#+I7+Ij 


We  can  now  use  the  identity  [5]: 


(100) 


(101) 


Q(a,p)  =  l-Q((3,a)  +  exp 


r_a2+p2 

v         2 


Io(a,P) 


(102) 


in  (101)  to  get 


Pr(R1<R2|l,Sj2)  =  Q 


V^p7        V2^ 


^2  +  4c+4j  V2  +  ^+^i 


exp 


Pc+Pj 


2VP^Pj 


1  +  ^c 


(103) 
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D.Pr(R1<R2|l,Sii,SJ2) 


This  probability  of  error  is  obtained  fromPrl  R,  <  R2|l,Sj   I  by  altering  the  signal 
in  branch  one  to  include  the  jammer  tone  S-  .  As  before,  a2c  is  replaced  with 

a2c  +  a2  j  +  2acaj  cos(0j  I .  We  must  then  integrate  twice  to  eliminate  the  dependence  on 

the  fading  parameters.  The  phase  angle  0  requires  a  third  integration.  However,  the 
integrations  cannot  be  evaluated  analytically.  The  triple  numerical  integration  is  time 
consuming.  As  will  be  explained  later,  an  effort  is  made  to  approximate  0  so  as  to  avoid 
the  third  integration. 

The  Prl  R,  <  R2  |l,Sj  )  for  a  non-fading  channel  is  given  in  [6]  as 


Pr(R,<R2|l,Sj2)  =  ^ 


1-Q 


"S 


ii  ~2    'A  ~2 
^  Y  a  n    Y  a  n  J 


+  Q 


'^ 


2      'A        2 
^  {  CT    n      Y  G    n  J 


(104) 


Using  (102),  we  get 


1-Q 


o2  o2 

a  c       a  ; 


2     M        2 
C7    n      Y  &    n 


=  Q 


a2c 


2      'A        2 
i^  V  °"    n      V  C7    n  ^ 


-exp 


<_a2c+aV 

V       2a2n    J 


^a   -aA 

UC  J 

^  cr2n  y 


(105) 


and  substituting  (105)  into  (104),  we  get 


PrR,<R2|l,S.     =Q 


7      ,V       2 
^  V  (7    n      \  G    n  J 


--exp 


^  o2  I      „2 

a  c  +a 


V       2a  n 


Mi, 


'a  -a ^ 

C  J 

V     CJ2n     J 


(106) 


Now,  replacing  a2c  with  a2c  +  a2j  +  2acaj  cos(0) ,  we  get 
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Pr(R1<R2|l,Sj,,Sj2,ac,aJ,0)  = 


a2j     fa2c +a2j +2acajCos(0) 


aV 


a2n 


f       -2 


exp 


a2c  +  2a2 j  +  2acaJ  cos(O) 
2a2n 


aj       a2c +a2j +2acaJcos(0) 


va- 


c\ 


(107) 
The  dependence  on  ac ,  aj ,  and  6  must  be  removed  by  multiplying  (107)  by  the  Ricean 

pdfs  for  the  signal  and  jammer  tones,  respectively,  and  by  the  pdf  for  0  and  then 
integrating  over  ac ,  a} ,  and  9.  As  previously  mentioned,  this  must  be  done  numerically 

and  can  be  a  computationally  intensive  operation.  In  any  event,  the  problem  at  hand  has 
been  "solved"  in  the  sense  expressed  in  [7]: 

In  many  cases  we  shall  consider  the  solution  found  if  we  can  express  it  in 
terms  of  one  or  more  integrals,  even  if  the  integration  can  only  be 
performed  by  expanding  in  series  (or  computing  numerically).  This  is 
called  'reducing  to  quadratures',  a  phrase  by  which  the  mathematician 
symbolically  washes  his  hands  of  the  remainder  of  the  task  of  finding  the 
solution. 

This  completes  the  analysis  of  the  receiver  performance  for  independent  multitone 
jamming.  The  final  form  of  the  equations  that  must  be  evaluated  are  summarized  in  the 
appendix,  equations  B.6  through  B.12.  Many  of  these  equations  are  also  used  for  the 
analysis  of  band  multitone  jamming,  which  is  addressed  in  the  next  chapter. 
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VI.  FH/MFSK  RECEIVER  PERFORMANCE  WITH  BAND  MULTITONE 

INTERFERENCE 


Band  multitone  jamming  implies  that  the  multiple,  equal  power  jamming  tones 
are  distributed  randomly  across  the  entire  frequency-hopped  bandwidth  with  a  specific 
number  of  jamming  tones  placed  in  a  frequency-hop  bin  given  that  a  particular  hop  bin  is 
jammed.  For  worst  case  band  multitone  jamming,  the  worst  probability  of  bit  error  is 
obtained  when  at  most  a  single  interfering  tone  exists  at  one  of  the  M  possible  orthogonal 
signaling  frequencies  in  a  specific  FH  band.  This  section  considers  the  general  case  of  a 
jammer  using  L  jamming  tones  per  hop  bin  with  1<L<M.  We  require  the  probability  that 
the  jammer  places  a  set  of  L  tones  within  the  signal  frequency-hop  bin.  This  probability  is 
derived  by  considering  the  following  numerical  example.  Suppose  that  we  have  a  total  of 
10  jammer  tones(q  =  10),  and  two  jamming  tones  per  jammed  bin  (L  =  2).  Suppose  that 
the  number  of  hop  bins  is  30(N=30).  The  jammer  tones  are  grouped  into  q/L  =  5  sets, 
each  set  composed  of  two  tones.  Now,  the  probability  that  one  of  the  five  sets  exists  in 
one  of  the  N  =  30  hop  bins  is  5/30.  In  the  general  case,  the  probability  that  the  hop  bin  is 
jammed  is 

Pr(bin  jammed)  =  ^—  =  -~~  .  (108) 

v  '       N        N-L 


The  probability  that  the  hop  bin  is  not  jammed  is 


Pr(binnotjammed)  =  l--V  (109) 

v  '  N-L 


From  the  total  probability  law,  we  find  that  the  overall  probability  of  error  is: 
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Pr(error)=  Pr(L,  L  >  l)  •  Ps  (L,  L  >  l)  +  Pr(L  =  0)  •  Ps  (L  =  0) 


N-L 


Ps(L,L>l)  +  (l-^-)-Ps(L  =  0) 


(110) 


The  probability  of  symbol  error  when  a  hop  bin  is  not  jammed,  PS(L  =  0) ,  is  the  same  as 
the  probability  of  symbol  error  for  noncoherent  MFSK  and  no  tone  interference 


Ps(L  =  0)<(M-l)-Pr(R,  <R2|1). 


(Ill) 


This  non-jamming  case  must  be  considered  in  order  to  account  for  every  possible 
jamming  scenario. 

The  probability  of  symbol  error  when  a  hop  bin  is  jammed  is  found  in  (38)  to  be 
(L>1): 


P*(L)^ 


+   1- 


(L-l)-Pr(R,  <R2|1,S.  ,S-  )  +  (M-L)-Pr(R,  <R2|l,SJi)] 
L-Pt(r,  <R2|l,Sj2)  +  (M-L-l)Pr(R,  <R2|l)j 


(112) 


Hence,  combining  (1 10),  (1 1 1)  and  (112),  we  obtain  the  probability  of  symbol  error 
when  (L>1)  as 
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Ps(L)  < 


1- 


N-L 
f 


M 


(L-l)-Pr(R,  <R2|l,SJi,SJi)+(M-L)-Pr(R1<R2|l,SJi) 
(l-^J-[L.Pi(R1<R2|l,SJi)+(M-L-0-Rr(Ri<Rili)' 


+ 


N-L 


(M-l) 


V 


2  +  4, 


exp 


\^ 


V    2  +  U) 


(113) 


The  individual  probabilities  appearing  in  (1 13)  are  identical  to  the  ones  found  in  the  last 
chapter,  and  are  given  in  (58),  (79),  (103),  and  (107).  The  difference  from  the 
independent  multitone  interference  is  that  L  is  fixed. 

This  completes  the  analysis  of  the  receiver  performance  for  band  multitone 
jamming.  The  final  form  of  the  equation  that  must  be  evaluated  is  given  in  the  appendix 
as  equation  B.  13.  In  the  next  chapter,  numerical  results  are  obtained  for  receiver 
performance  with  both  band  and  independent  multitone  jamming. 
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VII.  NUMERICAL  RESULTS 


The  performance  of  the  FH/MFSK  receiver  is  investigated  as  a  function  of  the 

P 
ratio  of  the  information  signal  power  to  the  total  interference  power,  — — .  Worst  case 

Jt 

performance  is  approximated  by  selecting  q  as  the  largest  integer  value  between  1  and  the 

Pjx 
ratio  — - .  For  non- worst-case  performance,  we  select  q=100. 

c 

The  non-fading,  Ricean,  or  Rayleigh  fading  channel  for  the  signal  and/or  the 
jammer  tones  is  modeled  by  appropriately  selecting  the  values  of  the  corresponding 

a2c  a2 

direct-to-diffuse  power  ratios,       2     and       [    .  A  large  ratio  value  implies  little  fading, 

2g  c  2a  j 

a  value  of  10  implies  Ricean  fading,  and  a  very  small  value  implies  Rayleigh  fading.  Due 

to  programming  necessities  and  in  order  to  avoid  divisions  by  zero,  we  use  a  ratio  equal 

to  10"9  to  simulate  Rayleigh  fading.  For  purposes  of  numerical  computation,  all  curves 
are  obtained  with  Eb/N0  =  13.35  dB  and  N=1000.  The  zero  thermal  noise  curve  is  also 
plotted  in  each  figure  to  validate  the  results  since  this  provides  an  upper  bound.  For  band 
multitone  interference  with  a  single  interfering  tone  per  hop  and  zero  thermal  noise,  the 
probability  of  bit  error  is  equal  to  [8] 


P   -± 


1    fPjT         1 
v—  -P. 

v  q       J 


2U 


(114) 
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while  for  independent  multitone  interference  with  zero  thermal  noise,  the  probability  of 
bit  error  for  FH/BFSK  is  [9] 


^      2Nl      2N-lJ    2   I  q        c 


+ 


q  r g-i y 


-i 


cos 


u 


I  q 


> 


(115) 


A.  BAND  MULTITONE  INTERFERENCE 

The  performance  of  FH/BFSK  with  band  multitone  interference  when  the 
information  signal  is  essentially  unaffected  by  channel  fading  (a2c  /  2a2c  =  1000)  with 

M=2  and  M=8,  respectively  is  shown  in  Figures  2  and  3.  L  successively  takes  the  values 

1  and  2  in  Figure  2  and  the  values  1,  4,  and  8  in  Figure  3.  Performance  is  computed  for 

Pc 
both  q  =  100  and  worst  case  q.  When  P3  >  Pc ,  or  equivalently  —  <  -20dB  and  q=100, 

JT 

the  worst  probability  of  bit  error  is  obtained  for  at  most  a  single  interfering  tone  at  one  of 
the  M  possible  orthogonal  signaling  frequencies  (L=l).  For  the  region  where  Pc  >  Pj ,  the 

worst  case  performance  is  obtained  when  all  signal  tones  are  jammed  (L=M).  Hence,  a 
jammer  whose  power  Pj  is  guaranteed  to  exceed  the  signal  power  will  do  more  harm  with 
just  one  tone  per  bin.  However,  a  jammer  with  Pj  guaranteed  to  be  less  then  the  power  of 
a  signal  tone,  will  do  more  harm  by  jamming  all  M  tones  of  a  bin.  For  worst  case  q,  the 
worst  performance  is  obtained  when  L=l  for  all  power  ratios. 

The  non-fading  signal  case  for  M=2  with  all  three  possibilities  of  jammer  tones 
fading  is  shown  in  Figure  4.  The  results  are  the  same  as  those  in  [1].  When  Pj  >  Pc  or 


equivalently 


< 


-20dB  and  q=100,  the  difference  in  performance  is  trivial  regardless 


Jt 


of  whether  or  not  the  jammer  tones  are  affected  by  the  fading  nature  of  the  channel,  even 
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when  the  jammer  tones  experience  Rayleigh  fading.  On  the  other  hand,  when  Pc  >  ?} ,  or 

Pc 

equivalently  — —  >  -20dB ,  the  interference  tones  that  experience  Rayleigh  fading  yield 

JT 

the  poorest  performance.  For  the  worst  case  q,  the  difference  is  trivial  for  the  entire  range 
of  signal-to-jammer  power  ratios.  Figure  4  is  identical  to  the  figure  we  get  if  we  use  the 
exact  result  equations.  In  fact,  the  individual  probabilities  of  bit  error  used  in  the  exact 
result,  as  they  are  described  in  [1],  are  identical  to  those  obtained  with  the  union  bound 
for  the  M=2  case.  So  M=2  will  not  be  treated  again  by  itself  since  it  can  convey  no  more 
information.  It  is  important  that  for  M=2  the  simpler  union  bound  equations  provide  the 
exact  result. 

The  non-fading  signal  case  for  M=8  is  shown  in  Figures  5  and  6  for  q=100  and 
worst  case  q  respectively.  The  curves  obtained  for  the  exact  result  and  the  union  bound 
are  identical  and  practically  indistinguishable  from  each  other.  The  previous  phenomenon 
with  regard  to  the  effect  of  fading  interference  tones  deduced  from  Figure  4  is  also 
obvious  here,  but  it  is  even  more  pronounced,  with  the  non-fading  and  Rayleigh  fading 
curves  differing  by  5dB. 

The  Ricean  fading  signal  case  for  q=100  and  worst  case  q  is  shown  in  Figures  7 

and  8  respectively,  where  the  exact  and  approximate  solution's  curves  obtained  for  M=8 

are  compared.  In  Figure  7,  we  notice  that  the  union  bound  curves  are  almost  identical  to 

Pc 
the  exact  result  ones  and  only  in  the  region  where  — —  >  -lOdBdo  they  appear  to  have  a 

JT 

constant  difference  of  1.3 8-1 0"4.  In  Figure  8,  the  union  bound  curves  behave  as  expected, 
being  practically  indistinguishable  from  the  exact  result  ones  for  Pc  <  P:  ,  while  in  the 

opposite  case  they  appear  to  have  a  slight  difference  with  a  maximum  value  equal  to 
1.38-10"4.  The  union  bound  curves  again  work  as  an  upper  bound  on  the  probability  of  bit 
error. 

Summarizing  the  results  up  to  this  point,  since  the  differences  detected  are  trivial 
ones,  for  all  practical  purposes  we  infer  that  for  M  even  larger  than  eight,    we  can  easily 
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resort  to  the  union  bound  equations  to  estimate  the  probability  of  bit  error  when  the  signal 
tones  are  either  not  affected  at  all  by  the  fading  channel  or  the  signal  tones  experience 
Ricean  fading. 

The  exact  result  and  the  union  bound  curves  for  a  Rayleigh  fading  channel 
affecting  the  signal  tones  and  the  three  possible  cases  of  fading  conditions  for  the  jammer 
tones  are  shown  in  Figure  9.  The  difference  between  the  two  groups  of  curves  attains  a 
maximum  value  equal  to  0.0038.  This  difference  between  the  exact  and  union  bound 
curves  is  the  greatest  noticed  up  to  this  point.  The  union  bound  curves  still  act  as  a  bound, 
being  higher  than  the  exact  result  ones,  but  the  accuracy  is  questionable.  This  is  examined 
further  in  the  following  figures  where  we  determine  if  increasing  the  number  of  symbols 
M  causes  a  decrease  in  the  probability  of  error  as  should  normally  happen. 

The  non-fading  and  Ricean  fading  signal  curves  in  conjunction  with  every 
possible  fading  situation  of  the  jammer  tones  is  displayed  in  Figures  10  through  14.  The 
value  of  M  is  sequentially  changed,  taking  the  values  2,  4,  and  8.  In  these  figures,  no 
anomaly  is  detected,  and  the  curves  behave  correctly  in  all  cases  with  the  probability  of 
error  decreasing  as  the  number  of  symbols  M  increases.  However,  the  next  three  figures, 
Figures  15,  16,  and  17,  display  the  Rayleigh  fading  signal  case  in  conjunction  with  all 
three  possible  fading  conditions  of  the  jammer  tones.  In  all  cases  the  results  obtained  are 
incorrect.  The  union  bound  still  acts  as  a  bound,  but  the  probability  of  error  erroneously 
increases  as  M  increases. 

Summarizing  the  union  bound  approach  for  band  multitone  interference,  we 
see  that  for  the  non-fading  and  the  Ricean  fading  signal  case,  the  union  bound  approach 
gives  very  accurate  results  which  are  perfect  in  the  first  case  and  deteriorate  some  in  the 
second  though  still  providing  acceptable  accuracy.  However,  when  the  signal  tones 
experience  Rayleigh  fading,  the  union  bound  approach  should  be  avoided  because  the 
accuracy  obtained  is  not  sufficient.  From  a  computational  point  of  view,  the  union  bound 
approach  is  far  simpler  than  the  exact  approach.  The  exact  approach  equations  include 
joint  pdfs  of  two  random  variables  expressed  with  double  integrals.  On  the  contrary,  the 
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union  bound  equations  are  expressed  as  the  sum  of  individual  pdfs  of  one  random 
variable  and  thus  they  are  expressed  in  terms  of  one  or  more  integrals. 


B.  INDEPENDENT  MULTITONE  INTERFERENCE 

In  independent  multitone  interference,  the  number  of  tones  existing  per  frequency 
hop  bin  can  vary  from  zero  to  M,  where  it  equates  to  the  number  of  signal  tones.  As 
mentioned  before,  the  calculation  of  the  probability  of  error  here  is  time  consuming  due 
to  the  inability  to  further  simplify  the  solution  which  contains  three  integrations. 
However,  one  integration  is  eliminated  for  the  case  where  either  the  signal  or  the  jammer 
tones  are  not  affected  by  the  fading  channel,  and  two  integrations  can  be  eliminated  for 
the  case  where  both  the  signal  and  jammer  tones  are  non-fading. 

The  M=2  case  for  a  non-fading  signal  and  various  jammer  fading  cases  when  we 
have  independent  multitone  interference  is  shown  in  Figure  18.  As  with  band  multitone 
interference,  the  results  obtained  for  M=2  are  identical  with  the  ones  given  by  the  exact 
solution.  This  is  because  the  exact  solution  equations,  described  in  [1],  again  coincide 
with  the  union  bound  for  M=2,  enabling  us  again  to  use  the  union  bound  with  perfect 
accuracy.  As  in  band  multitone  interference,  we  notice  the  performance  degradation  from 
the  non-fading  jammer  case  to  the  Rayleigh  fading  jammer  case  when  Pc  >  Pj . 

The  non-fading  signal  case  for  M=8  is  shown  in  Figure  19.  It  is  obvious  that  the 
union  bound  approach  is  a  valuable  tool  in  this  case  since  the  exact  result  approach 
requires  very  complicated  probabilistic  equations  with  the  complexity  further  increased  as 
M  increases.  However,  the  union  bound  approach  provides  us  with  a  concise  expression 
for  any  value  of  M.  The  only  means  of  comparing  the  curves  obtained  here  with  a  pre- 
existing result  is  to  get  a  feeling  from  the  comparable  ones  obtained  for  band  multitone 
interference  case  bearing  in  mind  the  result  obtained  in  [1],  that  the  probability  of  error 
for  the  band  multitone  interference  is  poorer  than  for  independent  multitone  interference. 
In  fact,  comparing  Figures  5  and  19,  we  see  that  they  have  the  same  shape,  differing  by  a 
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factor  of  roughly  eight.  Also,  comparing  Figures  4  and  18,  we  see  that  they  differ  by  a 
factor  of  roughly  equal  two.  The  difference  between  band  and  independent  multitone 
interference  is  roughly  equal  to  M. 

Figure  19  is  characterized  by  the  phenomenon  previously  described.  For 
Pc  >  Pj  a  channel  with  no  fading  of  the  interference  tones  provides  an  advantage  of  up  to 

5dB  over  one  with  Rayleigh  fading. 

The  Ricean  fading  signal  case  for  M=2  with  all  three  possibilities  of  fading  of  the 
jammer  tones  is  displayed  in  Figure  20.  As  for  non-fading  signal  tones,  when  q=100  we 
notice  a  difference  in  performance  for  the  non-fading  and  Rayleigh  fading  jammer  tones 
when  Pc  >  Pj .  That  difference,  however,  is  now  less  pronounced.  For  the  worst  case  q, 

the  difference  between  the  respective  curves  is  again  trivial.  Comparing  this  figure  with 
Figure  3,  we  again  notice  that  the  probability  of  error  is  less  by  a  factor  of  roughly  two. 

The  Ricean  fading  signal  case  for  M=8  with  all  three  possibilities  of  fading  of  the 
jammer  tones  is  displayed  in  Figure  21.  Comparing  this  figure  with  Figure  7,  we  notice 
that  the  probability  of  error  difference  is  eight.  This  confirms  the  previous  observation 
that  the  decrease  in  probability  of  error  when  we  shift  from  band  to  independent 
multitone  interference  is  roughly  equal  to  M. 

The  Rayleigh  fading  signal  case  for  M=2  and  M=8  with  non-fading  interference 
for  q  =  100  and  worst  case  q  is  shown  in  Figure  22.  The  obvious  problem  is  the  same  as 
the  one  observed  for  band  multitone  interference.  The  accuracy  achieved  is  not  sufficient 
and  predicts  an  increase  in  the  probability  of  bit  error  for  an  increase  in  the  number  of 
symbols  M.  Thus,  the  union  bound  approach  should  be  avoided  when  the  signal 
experiences  Rayleigh  fading. 

The  remaining  figures,  Figures  23,  24,  25,  and  26,  show  the  results  of 
approximating  6  rather  than  integrating  over  all  of  its  possible  values.  Remember  that  by 
0  we  denote  the  difference  between  the  phase  of  the  signal  and  the  phase  of  the  jammer 
tones.  This  is  modeled  as  a  uniformly  distributed  random  variable  over  [0,  2n].  As 
mentioned  before,  in  the  beginning  of  this  section,  for  the  case  of  more  than  one  jammer 


46 


tone  per  frequency  hop  bin,  we  must  perform  three  numerical  integrations  in  order  to 

calculate  the  probability  of  error:  one  for  the  amplitude  of  the  signal  tone  fading,  one  for 

the  amplitude  of  the  jammer  tone  fading,  and  one  for  the  angle  G.  To  this  must  be  added 

the  embedded  calculation  of  the  Marcum's  Q-function  that  is  also  defined  in  integral 

form.  It  would  be  very  convenient  to  find  a  way  to  eliminate  one  numerical  integration.  In 

fact,  the  integration  over  G  was  experimentally  found  to  be  the  most  time  consuming  one. 

After  trying  different  angle  values  to  get  an  accurate  approximation,  it  is  found  that  for 

7  .       .      . 

G  =  —  •  n  the  approximation  is  so  good  that  the  difference  between  the  exact  and 

approximated  traces  can  barely  be  detected.  The  exact  and  approximated  curves  for  a 
non-fading  signal  and  various  fading  of  the  jammer  tones  for  M=2  and  M=8,  respectively, 
are  shown  in  Figures  23  and  24.  The  exact  and  approximated  curves  for  a  Ricean  fading 
signal  and  various  fading  of  the  jammer  tones  for  M=2  and  M=8,  respectively,  are  shown 
in  Figures  25  and  26.  This  approximation  results  in  a  significant  reduction  in 
computational  complexity  and  time. 
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VIII.  CONCLUSION 

From  a  probabilistic  performance  point  of  view,  the  results  obtained  verify  those 
of  [1  ]  that  when  Pj  >  Pc  the  effect  of  fading  jammer  tones  is  very  small  when  the  signal 

tones  are  not  affected  by  fading.  When  Pc  >  Pj  performance  is  poorer  when  the 

interference  tones  experience  fading.  When  the  information  signal  experiences  fading,  the 
effect  of  fading  interference  tones  is  less. 

With  regard  to  the  union  bound  approach,  it  is  shown  to  be  very  accurate  for  non- 
fading  or  Ricean  fading  signal  tones,  while  for  Rayleigh  fading  signals  its  use  should  be 
avoided. 

The  application  of  the  union  bound  for  the  Ricean  fading  signal  in  the  case  of 
independent  multitone  interference  for  M=2  and  M=8,  that  was  not  treated  in  [1],  shows 
that  the  probability  of  error  advantage  when  the  jammer  uses  independent  multitone 
interference  instead  of  band  multitone  interference  is  equal  to  a  factor  of  M.  Hence,  when 
M  =  8,  independent  multitone  interference  gives  probabilities  of  error  eight  times  less 
than  does  band  multitone  interference. 

In  case  of  the  independent  multitone  interference,  the  calculations  involved  can  be 

7 
much  simplified  if  instead  of  integrating  over  all  possible  values  of  9  we  let  0  =  —  •  n , 

thus  eliminating  the  most  time-consuming  integration. 

In  summary,  the  union  bound  approach  taken  in  this  thesis  not  only  simplifies  the 
computational  complexity  involved  in  obtaining  numerical  results,  but  also  allows  us  to 
analytically  evaluate  situations  that  are  not  amenable  to  solution  using  an  exact  approach. 
The  principle  drawback  is  the  breakdown  of  the  union  bound  approach  when  the  signal 
experiences  Rayleigh  fading. 
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Figure  2:  Performance  of  a  FH/BFSK  noncoherent  receiver  with  essentially 
no  information  signal  fading  for  band  multitone  interference  with 
various  conditions  of  fading  of  the  multiple  interference  tones  for 
L=l  and  L=2. 
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Figure  3:  Performance  of  a  FH/8FSK  noncoherent  receiver  with  essentially 
no  information  signal  fading  for  band  multitone  interference  with 
various  conditions  of  fading  of  the  multiple  interference  tones  for 
L=1,L=4,  andL=8. 
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Figure  4:  Performance  of  a  FH/BFSK  noncoherent  receiver  with  essentially 
no  information  signal  fading  for  band  multitone  interference  with 
various  conditions  of  fading  of  the  multiple  interference  tones.  The 
results  are  identical  with  those  obtained  using  the  exact  result 
equations. 
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Figure  5 :  Performance  of  a  FH/8FSK  noncoherent  receiver  with  essentially 
no  information  signal  fading  for  band  multitone  interference  with 
various  conditions  of  fading  of  the  multiple  interference  tones,  q  = 
100,  using  both  the  exact  result  and  the  union  bound  equations.  The 
curves  obtained  are  identical  and  practically  indistinguishable. 


54 


-0- 

-o- 

-0- 


Pc/PjT(dB) 

Exact  result  -  jammer  tone  no  fading 
Exact  result  -  jammer  tone  Ricean  fading 
Exact  result  -  jammer  tone  Rayleigh  fading 
Union-bound  result  -  jammer  tone  no  fading 
Union-bound  result  -  jammer  tone  Ricean  fading 
Union-bound  result  -  jammer  tone  Rayleigh  fading 
Performance  for  zero  thermal  noise 


Figure  6:  Performane  of  a  FH/8FSK  noncoherent  receiver  with  essentially 
no  information  signal  fading  for  band  multitone  interference  with 
various  conditions  of  fading  of  the  multiple  interference  tones,  qwc 
using  both  the  exact  result  and  the  union  bound  equations.  The 
curves  obtained  are  identical  and  practically  indistinguishable. 
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Figure  7:  Performance  of  an  8FSK  noncoherent  receiver  with  Ricean  fading 
of  the  information  signal  for  band  multitone  interference  with 
various  conditions  of  fading  of  the  multiple  interference  tones, 
q  =  100,  using  both  the  exact  result  and  the  union  bound  equations. 
The  union  bound  curves  behave  as  expected,  being  a  little  higher 

Pc 

than  the  exact  solution  ones.  For  —  >  -lOdB  we  lose  a  little  bit  of 

PJT 
accuracy  as  the  union-bound  curves  have  the  same  shape  with  the 

exact  result  curves  but  they  have  a  constant  difference  of  1.36 10" 4. 
As  in  the  M=2  case,  there  is  no  essential  difference  between  no 
fading,  Ricean  fading  and  Rayleigh  fading  of  the  jamming  tones, 
for  Pj>Pc- 
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Figure  8:  Performance  of  an  8FSK  noncoherent  receiver  with  Ricean  fading 
of  the  information  signal  for  band  multitone  interference  with 
various  conditions  of  fading  of  the  multiple  interference  tones, 
qwc,  using  both  the  exact  result  and  the  union  bound  equations. 
The  union  bound  curves  behave  as  expected,  being  practically 
indistinguishable  from  the  exact  result  equations  for  Pc<Pjj  while 

for  Pc>Pjf  they  appear  a  little  bit  higher,  having  a  maximum 
difference  of  1.3410  4. 
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Figure  9:  Performance  of  an  8FSK  noncoherent  receiver  with  Ricean  fading 
of  the  information  signal  for  band  multitone  interference  with 
various  conditions  of  fading  of  the  multiple  interference  tones, 
q  =  100,  using  both  the  exact  result  and  the  union  bound  equations. 
Again  the  fading  behaviour  of  the  jammer  tone  is  irrelevant  to  the 
probability  of  error.  The  union  bound  curves  have  the  same  form  as 
the  ones  from  the  exact  result,  but  appear  to  display  the  biggest  error 
than  in  any  case  treated  up  to  is  point,  an  error  having  a  maximum 
equal  to  0.038 . 
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Figure  10:  Checking  the  Union  Bound  accuracy  for  essentially  no  fading  of  the 
information  signal  for  band  multitone  interference  with  no  fading  of 
the  jammer  tones,  for  both  q  =  100  and  qwc ,  for  M  =  2,  4,  and  8.  No 
anomaly  detected  for  the  Pb  using  the  union  bound  and  Pb  decreases 
as  M  increases. 
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Figure  1 1 :  Checking  the  Union  Bound  accuracy  for  essentially  no  fading  of  the 
information  signal  for  band  multitone  interference  with  Rayleigh 
fading  of  the  jammer  tones,  for  both  q  =  100  and  qwc ,  for  M  =  2,  4, 
and  8.  No  anomaly  detected  for  the  Pb  using  the  union  bound  and 
Pb  decreases  as  M  increases. 
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Figure  12:  Checking  the  Union  Bound  accuracy  for  Ricean  fading  of  the 

information  signal  for  band  multitone  interference  with  essentially 
no  fading  of  the  jammer  tones,  for  both  q  =  100  and  qwc ,  for  M  =  2. 
4,  and  8.  No  anomaly  detected  for  the  Pb  using  the  union  bound  and 
Pb  decreases  as  M  increases. 
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Figure  13:  Checking  the  Union  Bound  accuracy  for  Ricean  fading  of  the 
information  signal  for  band  multitone  interference  with  Ricean 
fading  of  the  jammer  tones,  for  both  q  =  100  and  qwc ,  for  M  =  2, 
4,  and  8.  No  anomaly  detected  for  the  Pb  using  the  union  bound 
and  Pb  decreases  as  M  increases. 
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Figure  14:  Checking  the  Union  Bound  accuracy  for  Ricean  fading  of  the 

information  signal  for  band  multitone  interference  with  Ray  lei  gh 
fading  of  the  jammer  tones,  for  both  q  =  100  and  qwc ,  for  M  =  2, 
4,  and  8.  No  anomaly  detected  for  the  Pb  using  the  union  bound 
and  Pb  decreases  as  M  increases. 
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Figure  15  :  Checking  the  Union  Bound  accuracy  for  Rayleigh  fading  of  the 
information  signal  for  band  multitone  interference  with 
essentially  no  fading  of  the  jammer  tones,  for  both  q  =  100  and 
qwc ,  for  M  =  2,  8,  and  16.  Anomaly  detected!!!  For  increasing  M 
the  probability  of  error  increases. 
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Figure  16  :  Checking  the  Union  Bound  accuracy  for  Rayleigh  fading  of  the 
information  signal  for  band  multitone  interference  with 
Ricean  fading  of  the  jammer  tones,  for  both  q  =  100  and  qwc , 
for  M  =  2,  8,  and  16.  Anomaly  detected! ! !  For  increasing  M  the 
probability  of  error  increases. 
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Figure  17:  Checking  the  Union  Bound  accuracy  for  Rayleigh  fading  of  the 
information  signal  for  band  multitone  interference  with 
Rayleigh  fading  of  the  jammer  tones,  for  both  q  =  100  and  qwc , 
for  M  =  2,  8,  and  16.  Anomaly  detected!!!  For  increasing  M  the 
probability  of  error  increases. 


66 


.Q 

CL 


10 


1=2 
Eb/No=13.35 


-30      -25      -20 


■15      -10       -5 
Pc/PjT(dB) 


0 


5 


10 


o    No  fading  jammer  tone 

*    Ricean  Fading  j  ammer  tone 

x    Rayleigh  Fading  jammer  tone 

+   No  fading  jammer  tone  worst  case  q 

~  Ricean  Fading  jammer  tone  worst  case  q 

-  -  Rayleigh  Fading  jammer  tone  worst  case  q 

. . .  Performance  for  zero  thermal  noise 

Figure  18:  Performance  of  a  FH/BFSK  noncoherent  receiver  with  essentially 
no  information  signal  fading  for  independent  multitone  interference 
with  various  conditions  of  fading  of  the  multiple  interference  tones. 
The  results  are  identical  with  the  ones  obtained  using  the  exact 
result  equations. 
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Figure  19:  Performance  of  a  FH/8FSK  non-coherent  receiver  with  essentially 
no  signal  fading  and  various  interference  fading  for  independent 
multitone  interference. 


68 


-10        -5 
Pc/PjT(dB) 


>   No  fading  jammer  tone 
:    Ricean  Fading  jammer  tone 
:    Rayleigh  Fading  jammer  tone 

-  No  fading  jammer  tone  worst  case  q 

-  Ricean  Fading  jammer  tone  worst  case  q 

-  Rayleigh  Fading  jammer  tone  worst  case  q 
. .  Performance  for  zero  thermal  noise 


Figure  20:  Performance  of  a  FH/BFSK  non-coherent  receiver  with  Ricean 
signal  fading  and  various  interference  fading  for  independent 
multitone  interference. 
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Figure  21 :  Performance  of  a  FH/8FSK  non-coherent  receiver  with  Ricean 
signal  fading  and  various  interference  fading  for  independent 
multitone  interference. 
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Figure  22:  Performance  of  a  FH/BFSK  and  FH/8FSK  receiver  for  Rayleigh 
fading  information  signal  and  non-fading  interference,  for 
independent  multitone  interference.  Anomaly  detected,  since 
for  an  increase  in  M  we  notice  an  increase  to  the  probability  of 
error. 
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Figure  23:  A  comparison  of  the  union  bound  curves  for  a  FH/BFSK 

noncoherent  receiver's  performance  with  essentially  no 

information  signal  fading  for  independent  multitone  interference 

with  either  no  fading  or  Rayleigh  fading  of  the  jammer  tones, 

with  q  =  1 00  and  worst  case  q,  obtained  by  integrating  9  and  also 

7 
by  approximating  0  with  —  n .  The  results  obtained  are  identical, 

so  the  approximation  can  effectively  be  used  simplifying  the 
calculations. 
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Figure  24:  A  comparison  of  the  union  bound  curves  for  a  FH/8FSK 

noncoherent  receiver's  performance  with  essentially  no 

information  signal  fading  for  independent  multitone  interference 

with  either  no  fading  or  Rayleigh  fading  of  the  jammer  tones, 

with  q  =  1 00  and  worst  case  q,  obtained  by  integrating  9  and  also 

7 
by  approximating  0  with  —n .  The  results  obtained  are  identical. 

so  the  approximation  can  effectively  be  used  simplifying  the 
calculations. 
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Figure  25:  A  comparison  of  the  union  bound  curves  for  a  FH/BFSK 

noncoherent  receiver's  performance  with  Ricean  information 

signal  fading  for  independent  multitone  interference  with  either 

no  fading  or  Rayleigh  fading  of  the  jammer  tones,  with  q  =  100 

and  worst  case  q,  obtained  by  integrating  8  and  also  by 

7 
approximating  0  with  —  n  .  The  results  obtained  are  identical,  so 

the  approximation  can  effectively  be  used  simplifying  the 
calculations. 
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Figure  26:  A  comparison  of  the  Union  Bound  Curves  for  a  FH/8FSK 

noncoherent  receiver's  performance  with  Ricean  information 

signal  fading  for  independent  multitone  interference  with  either 

no  fading  or  Rayleigh  fading  of  the  jammer  tones,  with  q  =  100 

and  worst  case  q,  obtained  by  integrating  0  and  also  by 

7 
approximating  0  with  — n  .  The  results  obtained  are  identical,  so 

the  approximation  can  effectively  be  used  simplifying  the 
calculations. 
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APPENDIX  .  MATHEMATICAL  IDENTITIES  AND  RELATIONS 


For  convenience,  all  pertinent  equations  are  repeated  here. 


A.  IDENTITIES 


1.  Integrals 


Mn(z)  =  exp[-^f].Jn(jz)  =  (-Jr.Jn(jz) 


(A.l) 


1     2n 
I0(z)  =  ^-  Jexp[±z-cos((j)-0)]d(l) 


(A.2) 


00 


m  +  n 


m!pnexp 


•  fx     2      -e_ax-J    (2(Wx1dx  = ±-? 

*  n^  K      '  m  +  n  +  1 
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J  ,-n 


m 


( *T\ 
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where  3nm(«)    is  the  Laguerre  polynomial: 
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\  av 


(A.3) 


00  , 

jx  •  exp[-  ax2  ]  •  I  n  (px)  •  J  „  (yx)dx  =  —  •  exp 


2        ..2 


4a 


•■£ 


(A.4) 
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2.    The  Marcum's  Q  function  and  related  integrals 


<    x2+a2^ 


r  x    i-a  f    x 

Q(a, P)  =  J  x  •  exp  -  — - —  •  I0 (ax)dx 


v         *      J 


(A.5) 


Q(a,0)  =  l 


(A.6) 


q(o,p)= 


exp 


(A.7) 


Q(a,p)  =  l-Q((3,a)  +  exp 


C       „2     ,    n2 


aJ_+P' 

^         2 


Io(a,P) 


(A.8) 
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G22 


CT2i  +CT22 
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a2, 


2  2     'A        2  2 

Q"i  +  G    2      \  <5    \  +G    2  J 

(A.9) 


B.  RELATIONS 


1.  Union  bound  on  the  probability  of  bit  error 

Ps  =  Pr(R,  <R2uR,<  R3u...uR,  <  RM|l) 

<Pr(R,  <R2|l)  +  Pr(R,  <  R3|l)+...+  Pr(R1  <  RM|1) 

<(M-l)Pr(R,  <R2|1) 
where  by  R,,  R2,...RM  we  express  the  random  variables  representing  the  outputs  of 
branches  1,2,...M,  respectively,  of  an  MFSK  receiver. 


(B.l) 
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2.    pdf  of  a  signal's  or  jammer's  amplitude  modeled  as  a  Ricean  random 
variable 


fAi(ai)  =  4--exp 


'    a2i+a2^  T  f  aVaV 


•u(ai)   0=c'j) 


(B.2) 


where  u(»)  is  the  unit  step  function,  I0(»)  is  the  modified  Bessel  function  of  the  first  kind 
and  order  zero,  a2c  is  the  power  of  the  direct  signal  component,  and  2a2c  is  the  power  of 
the  diffuse  signal  component  of  the  respective  tone. 


3.    pdf  of  a  signal's  or  jammer's  amplitude  modeled  as  a  Rayleigh  random 
variable 


fA;(a.)  =  ^r-exP 


^r}uW     0  =  c>i) 


a2, 


v    2c2 


(B.3) 


4.    Signal-to-noise  power  ratio  of  the  direct  signal  components 


a  i 

Pi  =-2-      (1  =  C>J) 
a  n 


(B.4) 


5.  Signal-to-noise  power  ratio  of  the  diffuse  signal  components 


2a2i 
£i=— —      (i  =  c,j) 

CJ    n 


(B.5) 
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6.  Pr(R,  <R2|l) 


Pr(R»<R2|l)  =  577-'exp 


2  +  ^( 


(B.6) 


7.  P^R^R^l.SjJ 


Pr(R,  <R2|l,Sj1)=         *   — -expl        Pc+Pj 


*  ii 


^j+4c+2       *\    4j+4c+2j     0Uj+4c+2J 


(B.7) 


8.  PrCR^RJ^S.) 


Pr(R1<R2|l,Sj2)  =  Q 
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(B.8) 


9.  P^R^RjfcS^Sj^.a^e) 


Pr(R1<R2|l,Sj,,Sj2,ac,aj,e)  = 


a2j       a2c +a2j+2acajcos(0)  |     l       (    a2c +2a2j +2acajCos(e) 


~2      '' 
CT    n 
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(B.9) 
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10.Pr(R1<R1|l,SJ1,Sj1) 


Pr(R1<R2|l,SJl,SJJ  =  — jf/J 

Z,'t    0  0  0 


2lt°0  00 
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(B.10) 


•fAe(ac)-fAj(aj)daedajde 


11.  Probability  of  symbol  error  for  L  jamming  tones  per  hop  bin,  1<L<M 


p*(L)4r 


+ 1- 


(L-l)-Ei{R1<Ra|I,Sjl.SJ|)+(M-L).Fr{R1<Rll|l>Sjb)]- 

L-Pt(r,  <R2|l,Sj2)  +  (M-L-l)-Pr(R,  <RM|l)] 


(B.ll) 


12.  Probability  of  symbol  error  for  independent  multitone  interference  for  L 
jamming  tones  per  hop  bin,  0<L<M. 


P.  =  infl - TJT^— r)l •  (M - 1) Pr(R,  <R2|l)  + 


N-M 

~L-1 


min(M-l.q) 

z 


L=l 
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M-N-i     J 
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where  the  probability  Ps(M)is  obtained  by  letting  L=M  in  PS(L).  This  yields 
Ps(M)  =  (M-l)-Pr(R1<R2|l,SJi,Sh). 
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13.  Probability  of  symbol  error  for  band  multifont  interference  for  L 
jamming  tones  per  hop  bin,  L>1 
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14.  Relation  between  probability  of  bit  error  and  probability  of  symbol  error 


A  yf   /  O 

P,  = •  Pc ,  where  M  is  the  number  of  symbols  (B.  1 4) 

M-l 


86 


INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center. 
8725  John  J.  Kingman  Rd.,  Ste  0944 
FT.  Belvoir,  VA  22060-6218 

2.  Dudley  Knox  Library 

Naval  Postgraduate  School 

411  DyerRd. 

Monterey,  CA  93943-5101 


Chairman,  Code  EC 

Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 
Monterey,  CA  93943-5121 

Prof.  R.  Clark  Robertson,  Code  EC/Re 

Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 
Monterey,  CA  93943-5121 


5.  Prof.  Tri.  T.  Ha,  Code  EC/Ha 

Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 

Monterey,  CA  93943-5121 

6.  Prof.  Roberto  Cristi,  Code  EC/Cx 

Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 

Monterey,  CA  93943-5121 


Prof.  Ralph  Hippenstiel,  Code  EC/Hi 

Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 
Monterey,  CA  93943-5121 


8.    Embassy  of  Greece 

Naval  Attache 

2228  Massachusetts  Ave.,  N.W. 

Washington,  DC  20008 


Lt.  JG.  George  Katsoulis. 
13PlaitouSt. 
Kallipolis,  Pireas,  18539 
GREECE 


87 


GUD'.EY  KNOX  LIBRARY 

"GRADUATE  SCHOOL 
MONI£Rt/CA    9SS43-5101 


